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Abstract 

We investigate the possibility of preparing left-handed materials in metallic mag- 
netic granular composites. Based on the effective medium approximation, we show 
that by incorporating metallic magnetic nanoparticles into an appropriate insulating 
matrix and controlling the directions of magnetization of metallic magnetic compo- 
nents and their volume fraction, it may be possible to prepare a composite medium 
of low eddy current loss which is left-handed for electromagnetic waves propagating in 
some special direction and polarization in a frequency region near the ferromagnetic 
resonance frequency. This composite may be easier to make on an industrial scale. 
In addition, its physical properties may be easily tuned by rotating the magnetization 
locally. 

PACS numbers: 73.20.Mf, 41.20.Jb, 42.70. Qs 



In classical electrodynamics, the response (typically frequency-dependent) of a ma- 
terial to electric and magnetic fields is characterized by two fundamental quantities, 
the permittivity e and the permeability /i. The permittivity relates the electric dis- 
placement field D to the electric field E through D = eE, and the permeability fi 
relates the magnetic field B and H hy B = iiH. If we do not take losses into account 
and treat e and // as real numbers, according to Maxwell's equations, electromagnetic 
waves can propagate through a material only if the index of refraction n, given by 
(e/i)^/^, is real. (Dissipation will add imaginary components to e and // and cause 
losses, but for a qualitative picture, one can ignore losses and treat e and n as real 
numbers. Also, strictly speaking, e and ^ are second-rank tensors, but they reduce 
to scalars for isotropic materials). In a medium withe and // both positive, the in- 
dex of refraction is real and electromagnetic waves can propagate. All our everyday 
transparent materials are such kind of media. In a medium where one of e and // is 
negative but the other is positive, the index of refraction is imaginary and electromag- 
netic waves cannot propagate. Metals and Earth's ionosphere are such kind of media. 
Metals and the ionosphere have free electrons that have a natural frequency-the plasma 
frequency-which is on the order of lOMHz in the ionosphere and falls at or above visible 
frequencies for most metals. At frequencies above the plasma frequency, e is positive 
and electromagnetic waves are transmitted. For lower frequencies, e becomes negative 
and the index of refraction is imaginary and consequently electromagnetic waves can- 
not propagate through. In fact, the electromagnetic response of metals in the visible 
and near ultraviolet regions is dominated by the negative epsilon concept[l]-[4]. 

Although all our everyday transparent materials have both positive e and positive /i, 
from the theoretical point of view, in a medium with e and both negative, the index 
of refraction is also positive and electromagnetic waves can also propagate through, 
moreover, if such media exist, the propagation of waves through them should give rise 
to several peculiar properties. This was first pointed out by Veselago over 30 thirty 
years ago- when no material with simultaneously negative e and // was known [5]. For 
example, the cross product of E and H for a plane wave in regular media gives the 
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direction of both propagation and energy flow, and the electric field E, the magnetic 

— * — * 

field H, and the wave vector k form a right-handed triplet of vectors. In contrast, 

— * — * 

in a medium with e and /i both negative, E x H ior a, plane wave still gives the 
direction of energy flow, but the wave itself (that is, the phase velocity) propagates 
in the opposite direction, i.e., wave vector k lies in the opposite direction of E x H 
for propagating waves. In this case, electric field E, magnetic field H, and wave 
vector k form a left-handed triplet of vectors. Such a medium is therefore termed 
left-handed medium[5]. In addition to this 'left-handed' characteristic, there are a 
number of other dramatically different propagation characteristics stemming from a 
simultaneous change of the signs of e and n, including reversal of both the Doppler 
shift and the Cerenkov radiation, anomalous refraction, and even reversal of radiation 
pressure to radiation tension. However, although these counterintuitive properties 
follow directly from Maxwell's equations-which still hold in these unusual materials, 
such type of left-handed materials have never been found in nature and these pecu- 
liar propagation properties have never been demonstrated experimentally. If such 
media can be prepared artificially, they will offer exciting opportunities to explore 
new physics and potential applications in the area of radiation-material interactions. 
Recently, interesting progress has been achieved in preparing a 'left-handed' material 
artificially. Following the suggestion of Pcndry[l], Smith and coworkers reported that 
a medium made up of an array of conducting nonmagnetic split ring resonators and 
continuous thin wires can have both an effective negative permittivity e and nega- 
tive permeability /i for electromagnetic waves propagating in some special direction 
and special polarization at microwave frequencies [6]. This is the first experimental 
realization of an artificial preparation of a left-handed material. Motivated by this 
progress, in this rapid communication, we propose to investigate the possibility of 
preparing left-handed materials in another type of systems-metallic magnetic granular 
composites. The idea is that, by incorporating metallic ferromagnetic nanoparticles 
into an appropriate insulating matrix and controlling the directions of magnetization 
of metallic magnetic particles and their volume fraction, it may be possible to achieve 
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a composite medium that has simultaneously negative e and negative // and low eddy 
current loss. This idea was based on the fact that on the one hand, the permittivity 
of metallic particles is automatically negative at frequencies less than the plasma fre- 
quency, and on the other hand, the effective permeability of ferromagnetic materials 
for electromagnetic waves propagating in some particular direction and polarization 
can be negative at frequency in the vicinity of the ferromagnetic resonance frequency 
cuq, which is usually in the frequency region of microwaves. So, if we can prepare 
a composite medium in which one component is both metallic and ferromagnetic and 
other component insulating, and we can control the directions of magnetization of 
metallic magnetic particles and their volume fraction, it may be possible to achieve 
a left-handed composite medium of low eddy current losses for electromagnetic waves 
propagating in some special direction and polarization. This composite may be easier 
to make on an industrial scale. In addition, its physical properties may be easily tuned 
by rotating the magnetization locally. 

To illustrate the above idea more clearly, in the following we present results of 
calculations based on the effective medium theory. Let us consider an idealized metallic 
magnetic granular composite consisting of two types of spherical particles, in which 
one type of particles are metallic ferromagnetic grains of radius Ri, the other type 
arc non- magnetic dielectric (insulating) grains of radius i?2- Each grain is assumed to 
be homogeneous. The directions of magnetization of all metallic magnetic grains are 
assumed to be in the same direction. In length scales much larger than the grain sizes, 
the composite can be considered as a homogeneous magnetic system. The permittivity 
and permeability of non-magnetic dielectric grains are both scalars, and will be denoted 
as ei and fii. The permittivity of metallic magnetic grains will be denoted as €2 and 
will be taken to have a Drude form: 62 = 1 — cUp/u^uj + i/r), where LUp is the plasma 
frequency of the metal and r is a relaxation time. Such a form of e is representative 
of a variety of metal composites [8] -[9]. The permeabihty of metallic magnetic grains 
are second-rank tensors and will be denoted as fi2, which can be derived from the 
Landau-Lifschitz equations[7]. Assuming that the directions of magnetization of all 
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magnetic grains are in the direction of the z-axis, jl2 will have the following form[7] : 



A*a -i/x' 

ill' jla 

1 



(1) 



where 



— * 

ujq — jHq is the ferromagnetic resonance frequency, Hq is the effective magnetic field 
in magnetic particles and may be a sum of the external magnetic field, the effective 
anisotropy field and the demagnetization field; cUm = jMq, where 7 is the gyromag- 
netic ratio, Mq is the saturation magnetization of magnetic particles; a is the magnetic 
damping coefficient; lu is the frequency of incident electromagnetic waves. We shall 
only consider incident electromagnetic waves propagating in the direction of the mag- 
netization. This is the most interesting case in the study of magneto-optical effects in 
ferromagnetic materials. We also assume that the grain sizes are much smaller com- 
pared with the characteristic wavelength A, and consequently, electromagnetic waves 
in the composite can be treated as propagating in a homogeneous magnetic system. 
According to Maxwell's equations, electromagnetic waves propagating in the direction 
of magnetization in a homogeneous magnetic material is either positive or negative 
transverse circularly polarized. If the composite can truly be treated as a homoge- 
neous magnetic system in the case of grain sizes much smaller than the characteristic 
wavelength, electric and magnetic fields in the composite should also be either positive 
or negative circularly polarized and can be expressed as : 

E{f,t) = E^^\ikz-0z-iujt 

where E^'' = x ^ iy, i?o^'' = x ^ iy, k = Real[A;e//] is the effective wave number, 
(3 = lm[keff] is the effective damping coefficient caused by the eddy current, kgff — 



k + ij3 is the effective propagation constant. In Eqs.(4)-(5) the signs of k and /? can 
both be positive or negative depending on the directions of the wave vector and the 
energy flow. For convenience we assume that the direction of energy flow is in the 
positive direction of the z axis, i.e., we assume /3 > in Eqs.(4)-(5), but the sign 
of k still can be positive or negative. In this case, if /c > 0, the phase velocity 
and energy flow are in the same directions, and from Maxwell's equation, one can 
see that the electric and magnetic field E and H and the wave vector k will form a 
right-handed triplet of vectors. This is the usual case for right-handed materials. In 
contrast, if < 0, the phase velocity and energy flow are in opposite directions, 
and E, H and k will form a left-handed triplet of vectors. This is just the peculiar 
case for left-handed materials. So, for incident waves of a given frequency cu, we can 
determine whether wave propagations in the composite is right-handed or left-handed 
through the relative sign changes of k and p. In the following, we shall determine 
the effective propagation constant keff = k + i(3 hj means of the effective medium 
approximation. The details of various kinds of effective medium approximations have 
been discussed in a series of references[8]-[12], here we only list the main points. First, 
if the composite can truly be considered as a homogeneous magnetic system in the case 
of grain sizes much smaller than the characteristic wavelength, then for waves (positive 
or negative circularly polarized) propagating through the composite in the direction of 
magnetization, their propagations can be described by an effective permittivity eg// 
and an effective permeability i^eff, which satisfy the following relations 

J D{f,uj)e'''^ff'df = Seff J E{f,uj)e'''^ff'df, (6) 
J B{r,u;)e'''^ff'dr = fXeff J h{f,u)e'''^ff'df, (7) 

where keff and u are related by keff — uj[eefffJ>effV^^- Although these relations are 
simple and in principle exact, it is very difficult to calculate the integrals in them 
because the fields in the composite are spatially varying in a random way. One 
therefore must resort to various types of approximations. The simplest approximation 
is the effective medium approximation. In this approximation, we calculate the fields 
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in each grain as if the grain were embedded in an effective medium of dielectric constant 
eg// and magnetic permeabihty ^Xeff■ Consider, for example, the ith grain. Under 
the embedding assumption, the electric and magnetic fields incident on the grain are 
the form of Eqs.(4)-(5): 

hinc = h^^^e'^^if'-'^^ (9) 

where E^'^ = x^iy and hl^^ = x^iy, corresponding to the positive(+) or negative(— ) 
circularly polarized waves. If the fields inside the grain can be found, then the inside 
fields can be used to calculate the integral over the grain volume 

li ^ [ Ei{f,uj)e'''^"'df, (10) 

Jvi 

Ji ^ [ hi{r,u)e'''^ff'df, (11) 

Jvi 

which is required to find the integral in Eqs.(6)-(7). For the positive or negative 

— * — * 

circularly polarized incident waves described by Eqs.(8)-(9), the integral /j and Jj can 
be written as 

li = {x^iy)Ii, (12) 
J- = {xTw)Ji, (13) 

where /j and Jj are scalars. If and J, can be found, then from Eqs.(6)-(7), the 
effective permittivity eg// and effective permeability /Xg// can be calculated by 

/lei/i + /2e2/2 



-eff 



/1/1 + /2/2 ' 



(14) 



where /i and /2 are the volume fractions of the two types of grains, //i is the perme- 
ability of non-magnetic dielectric grains, ix^^^ = ^a—fJ'' and /i^"-* = /ia+yu' (see Eqs.1-3) 
are the effective permeability of magnetic grains for positive and negative circularly 
polarized waves respectively. As for the calculation of li and Jj, we can follow the 



method of expanding interior and exterior fields in a multipole series and matcliing tlie 
boundary conditions [13]. After the coefficients of the multipole expansion of interior 
and exterior fields are obtained by matching the boundary conditions, /j and Ji can 
be found and subsequently be substituted into Eqs.(14)-(15). Since this method is 
standard, we shall not present the details. In the final results, Eqs.(14)-(15) reduce 
to one self-consistent equation: 

h Kffm^R^)akeffRi) - hiji{hR,)Q{KffR,) 

kiil)[{kiRi)i)i{keffRi) - KffiJi{hR,)ij'i{k^ffR,) ^ ^ 

kii^l{kiRi)Cl{keffRi)-keffMkiRi)Cl{keffRi)^ ' ^ ^ 

where Ri is the radius of the ith type of grains, and 

k^ = cu[ei/xi]i/2, (17) 
k2 = c^[62//?V/', (18) 

M^) = (19) 

Ci{x) = xh^\x), (20) 

ji{x) and hi{x) are the usual spherical Bessel and Hankel functions. Eq.(16) determine 
the effective product of {efi)eff, or equivalently ke/f, but not a single eg// and /ie//- It 
can describe the change of the phase of a plane wave across a slab of the composite, 
but it does not precisely describe wave propagations across a slab of the composite. 
This is due to the fact we make no attempt to rigorously solve the boundary-value 
problem for a slab of composite by matching the fields inside the slab and external 
fields outside the slab at the boundary. In fact, it is common in various types of 
effective medium theories that for u ^ the electromagnetic properties of a composite 
cannot in general be specified by a single eg// and i^eff- Since we can determine 
whether wave propagations through the composite is left-handed or right-handed by 
the calculation of the effective propagation constant keff, Eq.(16) is enough for the 
problems we are discussing. 

The numerical results for a metal volume fraction /2 of 0.3 obtained from Eq.(16) are 
summarized in Fig.l-Fig.2. Fig. 1(a) shows the frequency dependence of the real part of 



the effective permeability of magnetic grains for positive circularly polarized plane 
waves, Fig. 1(b) and (c) show the corresponding frequency dependences of the effective 
wave number k and the effective damping coefficient /3 in a composite consisting of 
metaUic magnetic grains and dielectric grains. Prom Eqs.(l)-(3), we can see that 
if the magnetic damping coefficient a is zero, Re[/x''+'*] will be negative in the whole 
frequency region of a; > a;o(the magnetic resonance frequency). From Fig. 1(a), we 
can see that if a is nonzero but small enough, there can still be a frequency region near 
LOo in which Re[//(+)] is negative. In this case, if the amplitude of the negative /i^^^ is 
large enough, k will be negative in this frequency region as was shown in Fig. 1(b), 
and hence the phase velocity and energy flow will be in the opposite directions in 
this frequency region, and E, H and k will form a left-handed triplet of vectors, 
i.e., the composite will be left-handed in this frequency region for positive circularly 
polarized plane waves. But if a is not small enough, Ref//^"^^] will be positive in 
the whole frequency region, or though Re[fi'^~^^] is negative in a frequency region near 
ujo, the amplitude of the negative Re[iJ,^^^] is not large enough, in this case k will be 
positive in the whole frequency region, as was shown in Fig. 1(b). In this case, the 
composite is right-handed for positive circularly polarized waves in the whole frequency 
region. The calculations also show that if the radius of metallic grains are small enough 
and the volume fraction of metal components is smaller than the threshold value of 
the insulator- metal transition, which is approximately 1/3 in our model, the losses 
caused by eddy current are very small and the composite is essentially an insulator. 
This can be seen from Fig.l(c), in which the damping coefficient /3 is very small 
compared with the amplitude of the wave number k, i.e., the eddy current losses are 
very small in the cases shown in Fig.l. If the volume fraction of metal components 
is larger than the threshold value, the composite will be essentially a metal, and the 
damping coefficient (3 will be much larger than the amplitude of wave number k{not 
shown in the figure). In Fig. 2 (a) we show the frequency dependence of the real part 
of the effective permeability of magnetic grains for negative circularly polarized 
waves, and in Fig. 2(b) we show the corresponding frequency dependence of the effective 
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wave number /c in a composite consisting of the metallic magnetic grains and dielectric 
grains. We can see that for negative circularly polarized waves, Re[/i^~)] is positive 
in the whole frequency region no matter how small a is, and correspondingly, k is 
positive in the whole frequency region, i.e., the composite is right-handed in the whole 
frequency region for negative circularly polarized waves no matter how small a is. 

In conclusion, we have discussed the possibility of preparing a left-handed material 
in metallic magnetic granular composites based on the effective medium approximation. 
Our model analysis shows that, by incorporating metallic magnetic nanoparticles into 
an appropriated insulating matrix and controlling the directions of magnetization of 
metallic magnetic components and their volume fraction, it may be possible to prepare 
a composite medium of low eddy current losses which is left-handed for electromagnetic 
waves propagating in some special direction and polarization in a frequency region 
near the magnetic resonance frequency. Further theoretical investigations by other 
approximations such as first-principle numerical calculations may be needed to further 
confirm the possibility shown in this paper. 
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Caption 



Fig.l. (a) The frequency dependencies of the effective permeabihty z/^"*"^ of magnetic 
grains and the corresponding frequency dependencies of (b) the effective wave number 
k and (c) the effective damping coefficient of the composite for positive circularly po- 
larized waves propagating in the direction of magnetization. (The plasma frequency cUp 
is usually in the visible or ultraviolet frequency region and the ferromagnetic resonance 
frequency ujq is usually in the microwave frequency region. For simplicity, hereafter 
wc will set uJo/ujp = 10~^. The other parameters are: cUm/i^o = 4.0, cUpR/c = 0.2, 
/2 = 0.3, a is shown in the figures). 

Fig. 2. (a) The frequency dependence of the effective permeabihty ii^"^ of magnetic 
grains and the corresponding frequency dependencies of (b) the effective wave number 
k and (c) the effective damping coefficient P of the composite for negative circularly 
polarized waves propagating in the direction of magnetization. (The parameters are: 
<^m/<^o = 4.0, ujpR/c — 0.2, /2 = 0.3, a is shown in the figures.) 
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